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ABSTRACT 

We have conducted a submillimetre mapping survey of faint, gravitationally lensed 
sources, where we have targeted twelve galaxy clusters and additionally the NTT 
Deep Field. The total area surveyed is 71.5arcmin^ in the image plane; correcting 
for gravitational lensing, the total area surveyed is 40 arcmin^ in the source plane for 
a typical source redshift z k, 2.5. In the deepest maps, an image plane depth of Icr 
r.m.s. ~ 0.8 mJy is reached. This survey is the largest survey to date to reach such 
depths. In total 59 sources were detected, including three multiply- imaged sources. 
The gravitational lensing makes it possible to detect sources with flux density below 
the blank field confusion limit. The lensing corrected fiuxes ranges from 0.11 mJy to 19 
mJy. After correcting for multiplicity there are 10 sources with fluxes < 2 mJy of which 
7 have sub-mJy fluxes, doubling the number of such sources known. Number counts 
are determined below the confusion limit. At 1 mJy the integrated number count is ~ 
10^ deg~^, and at 0.5 mJy it is ^ 2x 10* deg~^. Based on the number counts, at a source 
plan flux limit of 0.1 mJy, essentially all of the 850 /im background emission has been 
resolved. The dominant contribution (> 50 per cent) to the integrated background 
arises from sources with fluxes Ss.50 between 0.4 and 2.5 mJy, while the bright sources 
5*850 > 6mJy contribute only 10 per cent. 

Key words: Survey - submillimetre - galaxies: starburst - galaxies: high redshift - 
galaxies: evolution 



1 INTRODUCTION 

The first submillimetre mapping instrument SCUBA (Sub- 
millim etre Common User Bolometer Array; [Holland et al.l 
1 19991 ). which is mounted at the James Clerk Maxwell Tele- 
scope (JCMT) at Hawaii, allowed for observations of IR lu- 
minous gala xies at high re dshift. The first observations at 
850 /im (Sm all et al.l Il997h showed that these objects are 
much more common at earlier epochs. Subsequently, a num- 
ber of surveys have been undertaken to study this popu- 
lation of submm detected galaxies. The blank field surveys 
include observatio ns of t h e Hubble De ep Field Nort h (HDF- 
N) (. Hughes et al., 1 19981 : iBorvs et aL| 2003: Seric ant et all 
[iooa), the Hawaii Deep Fields (jBarger et al.lll999l'). C a nada- 
UK Deep SCUB A Survey (CUDSS) dEales et all I2OO0I : 
IWebb et al.1 |2003| ) the 8mJy survey (|Scott et all l2002h . 
Galactic regions (iBarnard et all |2004| ). the Groth strip 
dCoppin et"Sl l2005h . the SHADES survey (jCoppin et al.l 
I2OO6I ). and a re-analysis of several blank field surveys 
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(|Scott et al.ll2006l '). In particular CUDSS and SHADES have 
been successful in covering a large area of the sky. How- 
ever, the blank field surveys are limited by the confusion 
at 2mJy at 850 /im with the 15 m JCMT and hence do 
not probe the number counts of the fainter populatio n. The 
submillimetre extragala ctic background light (Fixs en et al.l 
1 19981 : iPuget et al.ll 19961 ) is, however, dom i nated by the pop- 
ulation around 1 mJy ( e-K- iBarger et al. 1 1 19991 : [ Blain et al.l 
ll999l : ICowie et al.ll200j ). To break the blank field confusion 
limit observing with SCUBA, gravitational lensing must be 
employed. The UK-SCUBA Lens Survey ((Small et al..il997l . 
I2OO2) targeted seven galaxy cluste r fields. Three of their 
fields were observed to larger depth ( Cowie et al. 20021). An - 



other lens survey was performed bv lChapman et al] ( 20021 ). 



however, this survey was relatively shallow. 

Submillimetre observations of objects at high redshifts, 
z > 1, benefit from the fact that the geometrical dimming of 
the light is cancelled by the negative fc-correction, resulting 
from the fact that the peak of the spectral energy distribu- 
tion (SED) is shifted toward the observing band. For a given 
luminosity, the observed submm flux is close to constant be- 



2 K.K. Knudsen, P.P. van der Werf, J. -P. Kneib 



tween redshift 1 and 8. Consequently, extragalactic submm 
observations primarily probe the high redshift universe. Fur- 
thermore, deeper surveys do not probe deeper into the uni- 
verse, but only sample lower luminosity galaxies. Galaxies 
in clusters at redshifts z < 1 are not expected to be seen 
with SCUBA, ex cept for sometimes the central cD galaxy 
lEdge et al.lll999l ') or an Active Galactic Nucleus (AGN). 

We here present the Leiden-SCUBA Lens Survey, in 
which we have targeted twelve galaxy clusters. This is the 
largest survey so far of gravitationally lensing clusters, and it 
is the first survey to substantially probe below the blank field 
confusion limit. This paper presents the observations, the 
analysis of the data, the resulting catalogue and the number 
counts. The analysis involves the mathematically rigo rous 
Mexican Hat wavelets algorithm (e.g. lCavon et al.l2000h and 
Monte Carlo simulations. In following papers, we will use 
the derived number counts as an observational constraint 
on models of the submm galaxy population, and we will 
present multiwavelength follow-up observations. 

In Section [2] we present the observations and the reduc- 
tion of the data. The source extraction is discussed in detail 
in Section |3] The issue of confusion is discussed in Section HI 
and the effect of gravitational lensing is discussed in Section 
[5l The resulting catalogue is presented in Section |6] Finally, 
in Section [7] we present the number counts for the survey. 
Throughout the paper we assume ilm ~ 0.3, 51a = 0.7 and 
Ho = 70km/s/Mpc. 



2 OBSERVATIONS AND REDUCTION 

We have obtained observations of a number of clusters of 
galaxies at 850 nm and 450 fim with SCUBA. In addition 
we ha ve obtained sim ilar observations of the NTT Deep 
Field (|Arnouts et al. I [1999.) . which was chosen due to the 
large, deep data set existing at optical and near-infrared 
wavelengths. In total our survey contains twelve fields of 
galaxy clusters and the one blank field covering an area of 
71.5arcmin^. The parameters for each field are listed in Ta- 
ble [1] 

SCUBA has two arrays of 37 and 91 bolometers opti- 
mized for 850 ^m respectively 450 /^m. A dichroic beamsplit- 
ter is used for simultaneous observations with both arrays. 
Both arrays have the bolometers arranged in a hexagonal 
pattern. Because SCUBA does not have a field rotator, the 
arrays appear as rotating on the sky. The field-of-view on 
the sky, which is approximately the same for both arrays, 
is roughly circular with a diameter of 2.3'. The observations 
were carried out in jiggle mode with a 64 point jiggle pat- 
tern, in order to fully sample the beam at both operating 
wavelengths. Subtraction of the strong sky background was 
done through 7.8 Hz chopping with the secondary mirror. 
Our observations were performed with a chop throw of 45" 
with the chopping position angle fixed in right ascension 
(RA). As a result the beam pattern has a central positive 
peak with negative sidelobes on each side, each with minus 
half the peak value, a pattern which can be used for the 
detection of at least the brighter sources. During the ob- 
servations the pointing was checked every hour by observ- 
ing bright blazars near the targeted fields. The noise level 
of the arrays was checked at least twice during an observ- 
ing shift, and the atmospheric opacity, r, was determined 



with JCMT at 850 fim and 450 /xm every two-three hours 
and supplemented with the T225GHz data from the neighbor- 
ing Caltech Submillimetre Observatory (CSO). Calibrators 
were observed every two to three hours. If available, primary 
calibrators, i.e. planets, preferably Uranus, were observed at 
least once during an observing shift. Our observations were 
supplemented with archival SCUBA data, hence the data 
set includes twelve cluster fields and the NTT Deep Field. 

The data were redu ced using the surf package 
(|jenness fc Lightfootlll998l ). First the chop of the secondary 
mirror was removed, i.e. the off-source measurements were 
subtracted from the on-source measurements. Then the 
varying responses of the bolometers were corrected by divid- 
ing with the array's flatfield. The extinction correction was 
performed based on the atmospheric opacity measured both 
with the JCMT and CSO. The and T4,50tJ.m was mea- 

sured a number of times during the night with the JCMT. As 
the atmospheric opacity may change on shorter timescales, 
the interpolated r-values may be somewhat inaccurate. At 
the CSO on the other hand, the opacity is measured ev- 
ery few minutes at 225 GHz. Using the linear relations be- 
tween r225GHz and rss pmn respectively T45o^m, deduced by 
[Archibald et al.l (|2000l ). it is possible to determine the atmo- 
spheric opacity at the time of the observations. The zenith 
opacity was for most of the time 0.12 < rgso/jm < 0.40. The 
data were inspected for bad or useless data. For each scan 
for each bolometer datapoints deviating by more than three 
sigma, based on the r.m.s. of the individual bolometer, were 
rejected from further analysis. This statistical exclusion is 
possible because there are no bright sources present in the 
data. Furthermore, the data were inspected by eye, and 
bolometers that were clearly more noisy than other bolome- 
ters were flagged and excluded from further analysis. The 
pointing of each scan was corrected using the pointing ob- 
servations taken just before and after the observations. The 
r.m.s. pointing error of the JCMT is typically 2". The corre- 
lated atmospheric fluctuations still present in the data were 
subtracted using the pixel-by-pixel median of the 25 least 
noisy bolometers. 

Data taken after January 2002, were affected by a pe- 
riodic noise of currently unknow n origin. This was flrst 
pointed o ut bv (jBor ys et al.^^2004 |), and has lat e r bee n dis- 
cussed bvlWebb et al. (2005), Sa wicki fc WebbI (I2005D . and 
ICoppin et al.l (120061. For these data, the power spectrum 
of the individual bolometers showed that some bolometers 
had a spike around l/16s. This spike, however, did not sys- 
tematically occur in the same bolometers or with the same 
strength. This effect was corrected for by performing a sky 
subtraction based on the bolometers that were not affected 
by this, and additionally, the affected bolometers were cor- 
rected through multiple-linear regression (T. Webb, private 
communication). As only a small fraction of the data for 
the survey was obtained after January 2002, this is rele- 
vant mostly for the MS 1054-03 data, where the data most 
affected was the Northern pointing. The correction for the 
l/16s spike brought down the noise in the affected data by 
~ 10 per cent. 

The scans were calibrated by multiplying by the ffux 
conversion factors (FCF) determined from the calibration 
maps. The FCFs are determined from the peak values (or 
the values corrected for extendedness) of the used calibra- 
tors. We estimate the uncertainty in the ffux calibration as 
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Table 1. The observed fields. 



NAME RA(J2000) Dee(J2000) z,, Csso <'oM, 

h m s o I II hours arcmin^ mJy/Cj, miy/C^ mjy/r2(, mJy/Qft 



C10016+16(*) 


00 


18 


33.2 


+16 


26 


17.8 


0.541 


7.73(5.46) 


4.5 


1.33 


2.00 


9.8 


16.5 


A478(*' 


04 


13 


25.3 


+10 


27 


54.3 


0.0881 


7.08 


4.3 


1.59 


2.05 


9.1 


14.5 


A496 


04 


33 


37.8 


-13 


15 


43.0 


0.0328 


10.4 


4.1 


1.08 


1.47 


11.6 


17.2 


A520 


04 


54 


07.0 


+02 


55 


12.0 


0.202 


19.6(18.0) 


4.3 


0.97 


1.26 


9.2 


14.5 


MS1054-03 


10 


56 


56.1 


-03 


36 


26.0 


0.826 


49.2 


14.4 


0.86 


1.49 


3.7 


10.2 


A1689(*) 


13 


11 


17.0 


-01 


20 


29.0 


0.181 


33.4(32.2) 


5.4 


0.70 


0.97 


4.4 


9.9 


RXJ1347.5-1145(*) 


13 


47 


30.5 


-11 


45 


09.0 


0.451 


10.5 


4.8 


2.04 


3.06 


7.2 


24.8 


MS1358+62 


13 


59 


50.6 


+62 


31 


05.1 


0.328 


4.80 


4.2 


1.39 


1.81 


7.6 


11.2 


A2204* 


16 


32 


46.9 


+05 


34 


33.0 


0.1523 


1.60 


4.0 


3.75 


5.20 


65.3 


87.2 


A2218 


16 


35 


54.2 


+66 


12 


37.0 


0.171 


42.3(35.6) 


7.7 


0.65 


1.06 


3.2 


16.0 


A2219(*) 


16 


40 


20.4 


+46 


42 


59.0 


0.225 


9.63 


4.6 


1.10 


1.54 


6.7 


11.8 


A2597 


23 


25 


19.8 


-12 


07 


26.4 


0.0852 


6.76 


4.1 


1.34 


1.77 


11.2 


17.8 


NTT Deep Field 


12 


05 


22.6 


-07 


44 


14.9 


-na- 


27.1 


5.0 


0.78 


0.97 


3.9 


6.3 



Parameters of the observed fields. The integration time tint is the total integration time, but without overheads (i.e., without the time 
needed for jiggle, chopping, etc). If the 450 /xm exposure time is different from the 850 /im its value is given in parentheses. The area, 
^850, given is the field covered after removing the noisy edge, a^^^^ is the lowest noise value in the whole field, cr^gf^td is the 
area-weighted noise level of the field, fli, is the beam. 

* Data from the JCMT archive. Supplemented with data from the JCMT archive. 



approximately 10 per cent at 850 pm. At 450 fj,m, the calibra- 
tion uncertainty is about 30 per cent, because of variations in 
the beam profile resulting from thermal deformations of the 
dish. This is in agreement with the canonical calibration un- 
certainties. The data were despiked by projecting the data 
on a grid and at each map pixel rejecting the associated 
bolometer pixels deviating by more than three sigma. Fi- 
nally, all the bolometers were weighted based on their mea- 
sured r.m.s. noise relative to one another and to the whole 
data set. The data were regridded with 1" pixels into a fi- 
nal map. The beam sizes are ~ 14.3" at 850 /im respectively 
~ 7.5" at 450 /im. The noisy edge was trimmed by remov- 
ing the outer 23 pixels; 23 pixels corresponds to one and a 
half beam. Unless otherwise mentioned the maps used in the 
analysis have been smoothed with a 5" full width half max- 
imum (FWHM) Gaussian to reduce high spatial frequency 
noise, resulting in beam sizes of 15.1" at 850 /im and 9" at 
450 fim. 



3 SOURCE EXTRACTION 

The sensitivity in the reduced maps is not uniform across 
the field. As we are working close to the signal-to-noise 
limit of these data it is crucial to understand the properties 
of the noise. Furthermore, the source extraction from maps 
like the SCUBA maps is a non-trivial task, which must be 
performed with a robust and well-un derstood method. In a 
previous paper (|Knudsen et aklboOSl ). we describe the ap- 
proach, which we adopt for this survey and first applied to 
the data for the field A2218. In summary, the noise is mea- 
sured across a map using Monte Carlo simulations, where 
the real data (i.e. the time streams) are substituted by the 
output from a random number generator with a Gaussian 
distribution and the same statistical properties as the real 
data. The simulated data is reduced following the same pro- 



cedure as the real data, creating empty maps, and the stan- 
dard deviation is measured for each pixel using about 500 
simulated maps. Furthermore, in order to remove the chop- 
ping pattern (caused by the motion of the secon dary mir- 
ror) from the beam we us e the CLEAN algo rithm (|Hogboml 
|l97i). For details, see iKnudsen et all l|2006l) . 



3.1 Point source detection 

Most sources at high redshifts have an angular extent on the 
sky of less than a few arcseconds. In the 850 fim beam such 
sources will appear as point sources. To do the point source 
extraction we choose to use Mexican Hat Wavelets (MHW). 
MHW is a mathematically rigorous tool for which the per- 
formance can be fully understood and quantified. MHW has 
proven to be a powerful source extraction technique in both 
SCUBA jiggle-maps and scan-maps (Barnard ot al. 20o3; 
iKnudsen etaLll2006l ). Wavelets are mathematical functions 
used for analysing according to scale. Isotropic wavelets have 
the advantage that no assumptions about the underlying 
field need to be made. The beam at 850 fim is well-described 
by a 2-dimensional Gaussian, which is best detected with 
the 'Mexican Hat' wavelets; the 'Mexican Hat' is the sec- 
ond derivative of a Gaussian. The software utilized are 
programmes written for the anticipated Planck Surveyo^ 
mission, but modified for application on SCUBA maps. 
For more details and the applicati o n of th e programmes 
on SCUBA maps see ICavon et al] (|2000D : iBarnard et al.l 
(|2004l ): lKnudsen et al.l (|2006l ). Here we summarize only the 
relevant details. 

The MHW source extraction is done in the following 
way. Point source candidates are selected at positions with 
wavelet coefficient values larger than a given number. For 
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each candidate, the wavelet function is compared to the the- 
oretical variations expected with the scale, as a further check 
on the source's shape. A value is calculated between the 
expected and the experimental results. If the is smaller 
than a given limit, i.e., the region surrounding the identified 
peak has the characteristics of a Gaussian point source of 
the correct dimensions, the point sour c e is in cluded in the 
extracted catalogue. In lKnudsen et ahl (|200 6!) we have per- 
formed controlled detection experiments to determine the 
optimal criteria for the MHW algorithm applied to SCUBA 
jiggle maps: the wavelet coefficient ^ 2 and x^ ^ 4. These 
we combine with the S/N ^ 3 criterion in the ffux map. 

In doing the source extraction, the MHW algorithm 
searches the maps for features at the scale of the beam, 
features with scales smaller or larger than the beam are 
not extracted. E.g., in the field RX J1347.5— 1145 the ex- 
tended emissi£n_£resumab^ due to the Sunyaev-Zel'dovich 



100 - 



Ply 

effect (|Komatsu et al.l [l999h is undetected by MHW point 
source detection. Also extended sources cannot be detected 
by changing the x^-hmit. This of course ensures that all 
sources detected are point sources. 

To check for undetected sources located closer than one 
beam to a detected source, the detected sources are sub- 
tracted from the map. If those are the only point sources 
present in the map, then MHW finds no sources in the resid- 
ual map. If any sources are significantly detected in the resid- 
ual map, those sources are subtracted from the original map 
and a new MHW detection is done in order to improve the 
result on the first detection. This is continued iteratively 
until the output converges. This approach does not resolve 

extended sources. 

As described in iKnudsen et"aLl (|2006l ). we performed de- 
tection experiments doing Monte Carlo simulations to anal- 
yse and determine the accuracy of the derived parameters 
(position and fiux). In the simulations point sources are 
added on to empty maps and recovered using the MHW 
algorithm. We perform this for all fields and find similar re- 
sults. The accuracies derived through this are included in 
the final errors quoted in Table [S] 



3.2 Completeness 

The completeness was determined through a similar set of 
simulations as those mentioned above, where point sources 
are added on to empty maps and recovered using the MHW 
algorithm with the same constraints as for the real data. 
This is done for each field, for a representative range of fluxes 
in steps of 1-2 mjy, repeated 4000 times for each flux level. 
The positions were chosen to be random with a uniform 
distribution. As the simulations are performed on the whole 
field, which has a non-uniform sensitivity, the results we 
find is an average across each field. In Figure [1] we plot the 
completeness for the individual fields. As is seen in Figure [TJ 
the completeness depends on the depth of the observations. 
For the deepest fields like A1689, A2218 and NTT Deep 
Field the observations are 80 per cent complete at a fiux 
level of ~ 3.5-4 mJy and 50 per cent at ~ 2.6-2.8 mJy. For 
the other, less deep fields the 80 per cent completeness is 
4.5-7.5 mJy and 50 per cent at 3.5-5.5 mJy. For A2204, which 
is a shallow field, the observations are 80 per cent complete 
at ISmJy and 50 per cent at 14mJy. 
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Figure 1. This figure shows the completeness as function of flux 
for all the surveyed fields. For a better overview the thirteen fields 
have been plotted in the three panels. The top panel shows the 
deepest of the fields, while the middle panel shows the medium- 
deep fields, and the lower panel shows two medium-deep fields, 
and shallower fields. Note that the x-axis is difli'ercnt for each 
panel. 



3.3 Spurious detections 

We have addressed the issue of spurious sources. Source de- 
tection was performed on inverted maps to check if any neg- 
ative sources were detected. We ignore negative sources sit- 
ting on a chop throw, which are a known artifacts from the 
original beam pattern. Furthermore, we ignore the deepest 
fields, where confusion plays a role. For those fields even 
a detected negative source can be a real structure in the 
background. Along the edges many negative sources were de- 
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tected. As this indicates that many positive spurious sources 
would be found along the edge, we have decided to trim 
away the edge at a width of one and a half beam, i.e. 23". 
In total we find five negative sources, which is in agreement 
with Gaussian statistics. Our catalogue thus may contain 5 
spurious sources. 

Additionally, we have performed source extraction from 
the Monte Carlo maps (see above) for all the fields. We find 
that in a hundred maps typically two sources were detected 
with S/N « 3, and no sources with S/N > 3.5. This demon- 
strates the reliability of MHW to pick out real sources with 
S/N > 3. With essentially no sources in the pure noise maps, 
the spurious sources detected in the inverted maps are prob- 
ably dAj^_to_tlM_stru£to background. 

ICoppin et al.l (|2006| ) used a Bayesian approach to fiag 
potential sp urious soures . In the method, which is described 
in detail in, Coppin et al.1 Hooi), an a priori probability dis- 
tribution is folded with the detection to calculate a pos- 
terior probability distribution, P{Sj), for each individual 
source. We adopt this approach and the prior is found 
determined from simulat ions of the submm sky assuming 
the number counts from (|Coppin et al.ll2006l ) and using the 
beam-pattern as known for the individual fields. We fiag the 
sources, which have a more than 5 per cent of their posterior 
probability distribution below OmJy, P{Sj < OmJy) > 5 
per cent. These sources will be marked in the catalogue ta- 
ble (Table[3|. The total number of sources is 11. The sources 
are not removed from the catalogue, as this is only a statis- 
tical approach and does not allow us to discriminate unam- 
bigiously between real and spurious sources as is evidenced 
in the field MS 1054-03. 



3.4 The 450 /im maps 

The atmospheric optical depth at 450 fim is five times larger 
than that at 850 nm. This makes it difficult to reliably cal- 
ibrate the 450 fim maps, as the data are much more sen- 
sitive to any variation in the sky opacity. As opposed to 
the 850 nm beam pattern, the 450 fim beam pattern is not 
well-described by a 2-dimensional Gaussian. In addition, the 
450 /im beam pattern is very sensitive to any deformation of 
the JCMT dish. Such deformations are normally the result 
of temperature variations. Consequently, the beam pattern 
changes during the night, with changing observing condi- 
tions. These effects complicate a reliable source extraction 
from the 450 fim. Furthermore, the 450 /tm beam is much 
narrower than the 850 /im beam, and more of the calibra- 
tion sources may appear extended. This has to be taken 
into account as well and might add to the uncertainty of the 
450 nm flux calibration. 

We have chosen not to perform deconvolution and 
MHW source extraction from the 450 fim maps to make a 
separate 450 /im point source catalogue. Instead we use the 
maps to determine the 450 fim flux at the 850 fim source po- 
sitions. That is done in the following way. The S/N map, 
the flux and the noisemap are used. In a circle with radius 
4" around the 850 fim source position, the maximum S/N 
value is found. If this value is equal to or greater than 3, 
the flux value in that pixel is adopted as the 5450 of the 
850 fjxa source, and the uncertainty, cts^^q , is 30 per cent of 
that value. If no pixels in that circle fuUfll the S/N criteria, 
an upper limit is given as three times the mean flux value 



corresponding to the three lowest noise values within the cir- 
cle. This has been done three times for all the 450 fim maps, 
namely when the maps have been smoothed with a Gaus- 
sian with FWHM of 5", 10" and 12.7". Smoothing with e.g. 
a 5" Gaussian reduces the high spatial frequency noise, while 
smoothing with a 12.7" Gaussian makes the 450 fim beam 
equivalent to the 850 fim beam. Furthermore, the smoothing 
reduces the effect that sources might be slightly extended 
in the original map due to intrinsic size and pointing er- 
rors. In a couple of cases the S/N criterion is met with one 
smoothing, but not quite in the other one(s). For the reasons 
described above, the 450 fim fluxes listed in Table |3] should 
used with caution. 



4 CONFUSION 

4.1 Confusion limit in blank fields 

In the background of deep SCUBA maps the instrumental 
noise and the confusion noise from a fainter submm pop- 
ulation are of a pproximately e qual magnitude. We adopt 
the formalism of lCondonI ll 19741 ') and use the same definition 



for the beam as iHogd (1 20011)7 ^beam = n{6 FW H M /2.35) , 
where 6fwhm is the width of the beam (SCUBA at 850 fim: 
Ofwhm = 15"). We adopt the rule of thumb that the num- 
ber of sources per beam should not exceed one source per 30 
beam s before the image is considered confused (e.g., iHogd 
I2OOII ). The confusion limit is the flux, Sconf, at which the 
number of sources in the map is one source per 30 beams. To 
estimate the confusion limit, the integrated number counts, 
A'^(> S), where A'^(> S) denotes the number density on the 
sky of sources brighter than S, are used. For blank fields 
a single power law suffices to describe the number counts, 
7V(> S) = NoS-°'. For N(> S) = l/30Qbeam, the confu- 
sion limit is Sconf ~ \/30flbeamNo. If we assume a = 2.0 
and A^o = 13000 deg"^ (based on e.g., iBarger et aLlll999l : 
iBorvs et al.ll2003t ). for SCUBA 850^01 blank field observa- 
tions, the confusion limit is ~ 2 mJy. For an average-sized, 
trimmed map in our survey of 5arcmin'^, this number of 
sources is Us — 4.7. This is close to the average number of 
sources per field in our survey. 

Confusi on in the ma ps affects the position and fiux de- 
termination. |Hog3 (| 200l l) have made general simulations ad- 
dressing the issue of the e rrors caused by confusion. Based 
on Figure 4 of his paper l|Hogell200lh . which gives the po- 
sition error as function of detected source density where no 
prior knowledge is used in the source detection, in the Eu- 
clidean case, a detected source density of one per 30 beams 
causes a median position error of 0.25 times the half width at 
half maximum (HWHM). At 850 fim this corresponds to an 
additional error in the position of 1.9". For crowded fields 
in our survey where the source density is larger than one 
per 30 beams, the error in the p osition caused by the con- 
fusion is 4-5". lEales et al. I l|2000l ) found in their simulations 
that in confusion-limited fields 10 per cent to 20 per cent of 
the detected sources would lie outside an error circle of 6". 
Furthermore, they found that the fluxes of the sources were 
boosted by a median factor of 1.44 , albeit with a large scat- 
ter. However, as argued bv iBlairil (|200ll ). confusion effects 
will only appear in SCUBA maps with detection limits of 
2 mJy or less at 850 ^m; hence most of our data are relatively 
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unaffected by flux boosting, though flux boosting is expected 
to play a role in the deepest maps, i.e. A1689, A2218, NDF 
and possibly A520. For these field s we adopt the Bayesian 
approach as used bv lCoppin et al. I (|2005. . .2006.), and" as al- 
ready mentioned in Subsect. 13.31 to estimate how large the 
flux boosting might be. The estimated de-boosted flux is 
listed in the catalogue Table |3l The deboosted fluxes are 
about 1 mjy fainter, though the actual number is connected 
with the S/N of the detection. 

4.2 Confusion limit in lensed fields 

For the cluster flelds the confusion limit is affected by the 
gravitational lensing. The gravitational lensing magnifles the 
region seen behind the cluster, hence the source plane is 
smaller than the image plane. The number of beams is con- 
served between the image plane and the source plane, i.e. 
the size of the beam scales with the magnification. This is 
why it is at all possible to observe the fainter sources, which 
have a higher surface density than the brighter sources. 

The number counts in the lensed case can be writ- 
ten as Niena = {No / ^i) {S / fi)~°' = Nbiank^"~^ , where 
fjL is the gravitational lensing magnification. The confu- 
sion limit in the lensed c ase can thus be written as 
Sconf = y^30QbeamNofi^~°' ■ As the Icusiug magnification 
varies across the field. We use the average magnification for 
a field as estimated by the ratio of the area in the image 
plane and the area in the source plane. This simple estimate 
gives an average across the field, which in some cases means 
that the estimated lensed confusion limit does not refiect the 
confusion limit of the highly magnified region close to the 
caustics. For the most extreme case in our sample, A1689, 
where the source plane area surveyed is 20 times smaller 
than the SCUBA field of view, the confusion limit is reduced 
by a factor 4.5, i.e., Sconf ~ 0.44 mJy. The confusion lim- 
its for the cluster fields based on this simple calculation are 
given in Table [2l In the simplified estimate of the confusion 
limit presented here we assumed that the number counts are 
described by a single power law. There are good indications 
that the number counts are described by a double power law 
or another function with a (gradual) turn-over (see Sect.[71). 
Including this in such a calculation will work in a favourable 
direction and the confusion limit in the source plane will be 
lower. 



5 GRAVITATIONAL LENSING 

To quantify the gravitationa l lensing effect we use 
LENSTOOL (|Kneib et al.l 119931 ). The gravitational poten- 
tial of each galaxy cluster is mapped in a mass model, which 
describes the distribution of the overall potential of the clus- 
ter and to some extent of the individual galaxies. For clus- 
ters where the cluster lensing is a less strong effect only the 
brightest cluster galaxies are considered in the mass model 
in addition to the global cluster potential. For clusters where 
the cluster lensing is a dominant effect many galaxies have 
been included to map the substructure of the total poten- 
tial, as individual galaxies might cause extra lensing of the 
background sources. The lensing correction is done for the 
individual submm sources and for the sensitivity maps. The 
latter are needed in the calculation of the number counts. 
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Figure 2. Histogram showing the distribution of the magnifi- 
cation factors. The histogram is binned with = 1, except 
between 1 and 2, where the bin-size is 0.5. 
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Figure 3. The area as function of magnification for the individ- 
ual fields. This has been split into two figures for clarity. The top 
panel shows the fields with a relatively less strong lensing effect, 
while the bottom panel shows those with a stronger lensing effect. 
We have assumed a source redshift of 2: = 2.5. Piacing the source 
plane at a different redshift, z > 1, would not make a significant 
difference. 

As the redshift is not yet known for a majority of the ob- 
jects, we assume z — 2.5 for the objects wit h unknown red- 
shift based on the redshift distribution from lChapman et al.l 
(2003, 2005). Likewise, the sensitivity maps, which give us 
the observational sensitivity in the image plane, are traced 
to a source plane at z = 2.5. The actual redshift distri- 
bution of the faint SMGs is not known and it is also not 
known whether it follows that of the brighter SMGs as de- 
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du ced by Chapman et alj l|2005h . Based on a stacking analy- 
sis, IWang et al.l l|2006h suggest that the redshift distribution 
of faint SMGs pe a ks at redshifts of one or less. However, 
m iKnudsen etal] l|200d ) submm stacking results of high 
redshift red galaxies show that half of the EBL light pro- 
duced by at the faint end originates from red galaxies in the 
redshift interval 1-4. Of the five /gso < 1 mJy SMGs with 
reliable identification and spectroscopic redshifts, the red- 
shifts are z ^ 1.0 2.5, 2.5, 2.6 and 2 9 l|Kneib et al.ll2004 
iBorvs et al ] |2004l : iKnudsen et aLll2006l . and see appendix for 
A1689), showing no evidence fo r a radically different d istri- 
bution than the one deduced bv lChapman et al.l l|2005l ). We 
note that not knowing the exact redshifts of the SMGs will 
introduce some uncertainty in the lensing correction, how- 
ever, the magnification correction is only weakly dependent 
on redshift at z > 1, where we expect most of the sources 
to be. 

The magnification factors of the individual sources 
range from 1.1 to 23. We have plotted a histogram of the 
magnification factors in Fig. [2] About 40 per cent of the 
sources are magnified by factors fi > 2, while 20 per cent are 
magnified by 1.5 < fj, < 2, and 40 per cent have relatively 
low magnification factors 1 < /i < 1.5. We have plotted the 
area as function of magnification factor for the individual 
fields in Figure [S] and the area as function of source plane 
sensitivity for the whole survey Figure [l] In Figure [S] the 
two most extreme cases are A1689, where the area surveyed 
is a factor 20 smaller in the source plane than in the im- 
age plane, and MS1054-03, where the average over the large 
angular area surveyed of the cluster dilutes the strong lens- 
ing effect caused by the core of the cluster. The total area 
observed by our survey in these 13 fields is 71.5 arcmin^. 
When taking into account the gravitational magnification 
the area of the cluster fields is reduced to 35 arcmin^ in the 
source plane. The area of the individual fields as well as the 
sensitivity in the source plane is listed Table [5] For com- 
parison the seven cluster fields from the UK-SCUBA Lens 
Survey are 40 arc min^ in the image plane and 15 arcmin^ in 
the source plane (ISmail et a l.' 2002). Likewise, in the deep 
though small survey by Jcowie ct al. 2002;), the area in the 
image plane is 18 arcmin^; assuming a reasonable amplifi- 
cation this corresponds to 6 arcmin^ in the source plane. 

The uncertainties of the corrected fiuxes and positions 
introduced by the lensing are in most cases small. The mag- 
nification is generally a monotonic function of the redshift 
(except in the very central part of a strong lensing cluster), 
but for source redshifts twice larger than the lens redshift, 
the amplification is only weakly increasing with redshift. As 
essentially all the SCUBA sources are expected to be at 
z >1, redshift dependence in the lensing correction is only a 
minor effect. However, the position of the source relative to 
the cluster centre and cluster members can play an impor- 
tant role, as the magnification can vary from ~ 2 to ~ 20 for 
typical background sources. The uncertainty in the magnifi- 
cation (assuming a known redshift) is directly linked to the 
uncertainty in the mass model. The closer the object is to 
a critical line, the higher its magnification will be and the 
larger the uncertainty on the magnification will be. At most 
four of the SCUBA sources lie relatively close to critical 
lines and therefore their magnification factors are subject to 
larger uncertainties. However, on an ensemble basis, when 




1 10 
S(850 /xm) [mJy] 

Figure 4. The area as function of Icr sensitivity for the whole 
survey. The dashed line represents the area in the image plane, 
and the solid line represents the area in a source plane at redshift 
z = 2.5. The actual area surveyed in the source plane is a factor 
two smaller than what we observe in the field of view of SCUBA 
at 850 iim. For the twelve cluster fields the total area of useful 
data in the image plane is 66.5 arcmin^ and in a source plane at 
redshift z = 2.5 is 35 arcmin^. 



e.g. deriving the counts in terms of unlensed flux, the error 
in magnification is compensated by the error in lensed area, 
thus the change in the unlensed counts due to uncertainty 
in the mass model is of the second order. 

We estimate the uncertainty of the magnification fac- 
tors of the individual sources through a Monte Carlo sim- 
ulation. The magnification is determined at 1600 positions 
with a normal distribution within the la error circle cen- 
tered on the MHW position. In Table [3] we give the median 
magnification from the MC simulations together with the 
68 per cent deviation of the magnifications determined at 
the MC positions. For very large magnification factors, typ- 
ically > 10, such as seen for the multiple-imaged galaxies in 
A1689, which are close to the caustics, the MC magnification 
factor distribution has both a large skewness and kurtosis. 
While we have made an estimate of the strength of possible 
fiux-boosting, the results of these MC simulations show that 
the uncertainties on the magnification factor often exceeds 
the flux-boosting. 

Identification of multiply-imaged galaxies in the sam- 
ple is important, as a repeated counting of the same source 
will affect the number counts. In total three multiply-imaged 
sources are found in the fields of the strongly lensing clus- 
ters A1689 and A2218. The galaxy in A2218 is triple- 
imaged with a total magnification factor of 45 and has been 
studied extensively Jkneib et all |2004| : ISheth et al] |2004| : 
iKneib et aLll2005l : [Garrett et al.ll2005l ). The two other galax- 
ies are present in A1689, one is triple-imaged and the other 
is a quintuple-imaged galaxy, both with spectroscopic red- 
shift ~ 2.5. Their identification will be discussed in a future 
paper. 



6 THE CATALOGUE 

In the twelve cluster fields we detect 54 sources and in the 
NTT Deep Field we detect 5 sources. The sources have been 
named according to their detection (SMM) and their J2000 
coordinates. The catalogue of the extracted point sources is 
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Table 3. Catalogue of source positions, submm fluxes and uncertainties. S is the flux, S/N is the signal-to-noise of the detection in flux 
units, the it's give the uncertainties in the flux and position. The uncertainties do not take into account possible additional uncertainties 
due to confusion. For the position, the additional confusion-uncertainty is ~ 1.9" (as described in Section|4)l. The * after the source name 
indicates that the source has a P{Sj < Omjy) > 5% (as described in subsection 13.31 1. fi is the lensing magnification, while ^mc is the 
lensing magnification from Monte Carlo simulations (see Section [5J. 585o(deboost) is the deboosted flux (see subsection 14. H . 



NAME 


O-pos Ss50 S/N O"S850 


5450 S/N4i,o CTs^gg 


z 






S'85o(deboost) 




" mjy mjy 


mjy mjy 








mJy 



CW016+16 



SMM J001828.9+162617* 


4.0 


6.5 


3.2 


2.2 


< 34.2 






2.5 


1 


2 


1.2 ± 0.02 


SMM J001829.4+162653* 


4.0 


5.8 


3.2 


2.0 


< 77.9 






2.5 


1 


2 


1.2 ± 0.02 


SMM J001834.2+162517 


4.0 


7.0 


3.9 


2.4 


< 56.1 






2.5 


1 


3 


1.3 ±0.03 


SMM J001835.H-162559* 


4.0 


5.3 


3.1 


1.8 


< 34.2 






2.5 


2 





2 0+"-^ 


AJ^78 
























SMM J041322.9+102806* 


3.8 


6.8 


3.3 


2.5 


< 35.3 






2.5 


1 


2 


1.2 ± 0.02 


SMM J041323.4+102657 


3.1 


7.9 


4.2 


2.1 


< 53.9 






2.5 


1 


3 


1.3 ±0.02 


SMM J041327.2+102743 


2.3 


25.0 


14.4 


2.8 


55.4 


5.3 


16.6 


2.837 


1 


3 


1.3 ±0.03 


SMM J041328.7+102805* 


3.8 


9.0 


3.8 


3.3 


< 32.8 






2.5 


1 


3 


1.3 ±0.02 


AJ^96 
























SMM J043334.7- 131526* 


4.1 


4.7 


3.1 


1.7 


< 43.6 






2.5 


1 


4 


1.4 ±0.03 


SMM J043335.4- 131454* 


4.1 


5.3 


3.3 


1.9 


< 54.7 






2.5 


1 


4 


1.4 ±0.02 


SMM J043336.5-131547 


3.2 


4.8 


4.0 


1.7 


51.1 


3.0 


15.3 


0.03 








SMM J043337.4-131558 


4.1 


4.7 


3.8 


1.7 


< 42.9 






0.03 








SMM J043337.6- 131627 


3.0 


9.0 


5.5 


1.6 


< 40.3 






2.5 


1 


5 


1.5 ± 0.03 


SMM J043337.8- 131541 


3.0 


7.9 


5.7 


1.4 


< 39.4 






0.03 








SMM J043338.9-131444 


4.1 


4.0 


3.1 


1.4 


< 67.3 






2.5 


1 


4 


1.5 ± 0.04 


SMM J043339.4-131637 


4.1 


4.8 


3.6 


1.7 


< 52.1 






2.5 


1 


4 


1.4 ±0.02 


SMMJ043340.1-131533 


3.0 


6.4 


5.1 


1.1 


< 36.9 






2.5 


1 


5 


1.5 ± 0.04 



A520 



SMMJ045403.1+025547 


3.3 


4.7 


4.1 


1.1 


< 31.9 






2.5 


1.5 


1.5 ± 0.03 


3.7± 1.3 


SMM J045406.2+025410* 


4.2 


3.9 


3.1 


1.4 


< 32.6 






2.5 


5.5 


5.5 ±0.4 


3.4 ± 1.3 


SMM J045406.7+025435* 


4.2 


4.3 


3.1 


1.5 


< 37.0 






2.5 


4.5 




2.3 ± 1.6 


SMM J045409.7+025510 


3.3 


6.0 


4.4 


1.4 


29.0 


2.7 


8.7 


2.5 


1.7 


1.7 ±0.05 


4.7 ± 1.5 


MS 1054- 03 
























SMM J105649.3-033606 


3.3 


5.0 


3.6 


1.1 


< 17.6 






2.5 


1.1 


1.07 ±0.003 




SMM J105655.8-033610 


3.3 


3.9 


3.8 


0.9 


25.6 


3.4 


7.7 


2.5 


1.1 


1.11 ± 0.007 




SMM J105656.3-033635 


3.3 


3.9 


3.7 


0.9 


< 21.2 






2.5 


1.1 


1.18 ±0.01 




SMM J105657.0-033612 


2.8 


4.9 


4.8 


0.9 


61.7 


3.6 


18.5 


2.5 


1.1 


1.12 ± 0.008 




SMM J105700.3-033513 


3.9 


3.5 


3.2 


1.1 


< 20.2 






2.5 


1.1 


1.05 ± 0.003 




SMM J105700.3-033544 


3.3 


4.4 


3.5 


1.0 


28.1 


3.6 


8.4 


2.5 


1.1 


1.08 ±0.004 




SMM J105701. 8-033827 


3.3 


4.7 


3.5 


1.1 


< 27.9 






2.5 


1.3 


1.3 ±0.03 




SMM J105702.2-033604* 


3.9 


4.4 


3.0 


1.4 


< 21.8 






2.423 


1.1 


1.11 ± 0.007 




SMM J105703.7-033730 


3.9 


4.2 


3.3 


1.4 


62.7 


5.7 


18.8 


2.5 


1.6 


1.6 ±0.05 





Continues on the next page 



given in Table [3] After correcting for lensing multiplicity, 
we have detected 15 sources below the blank field confu- 
sion limit. Of these, 7 have flux densities < 1 mJy, which 
doubles the number of known sub -mj y sources (compar e 
ICowie et al1l2002l . [sinafl et al.ll2002l and lBorvs et al.ll2004l ). 
A description of the individual fields can be found in Ap- 
pendix |A] along with the final maps. 



7 NUMBER COUNTS 

T.l Determining the number counts from maps of 
non-uniform sensitivity 

The notation 7V(> 5*) is typically used for cumulative num- 
ber counts: the number of sources per unit solid angle 



brighter than a flux limit S. Calculating the cumulative 
number counts by counting the number of sources with > S 
must be done on a map with uniform sensitivity. SCUBA 
maps, however, do not have uniform sensitivity. The prob- 
lem of determining number counts for maps of non-uniform 
sensitivity has previo usly been discussed for several of the 
blank fleld surveys JWebb et al.1 [2OO3I : iBorvs et all I2OO3I : 
ICoppin et alll2006l : [Scott et al.ll200d ). The presence of grav- 
itational lensing results in even larger non-uniformity com- 
pared to some of the large blank fleld surveys and com- 
plicates any completeness corrections. We here present the 
number counts as deduced using two different approaches. In 
both cases cluster members like cD galaxies are excluded and 
also the sources, which we have marked as potentially spu- 
rious using the scheme from iCoppin et al., (,2005. ) (see sub- 
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NAME 


ff 


5850 
mJy 


S/N 


mJy 


S450 
mJy 


5/Af450 


""5450 

mJy 


z 




MMC 


S85o(deboost) 


A 1689 
























SMM J131 125.7—012117 


3.3 


5.0 


3.9 


1.6 


< 66.4 






2.5 


3.9 


3.9 ±0.4 


3.7 ± 1.5 


SMM J131128. 6— 012036 


4.3 


2.6 


3.4 


0.8 








2.5 


23.6 


17 6+?S1 

^' -"-10.6 


1.9 ± 0.9 


SMM J131128. 8-012138 


4.3 


3.6 


3.3 


1.2 


< 32.9 






2.5 


5.8 




2.4 ± 1.1 


Oi.vi.ivi o ±i) 1- . J- — yj A.^\j'-tiy 


2.8 


4.7 


6 


0.8 


21.4 


4.4 


6.4 


2.5 


21.6 


IS 7+19-1 


4.3 ± 0.8 


^A/TA/T T1 'i^ 1 9Q S 01 90*^7 




9 ^ 




U.o 


^ i 1 .0 






Ci. 


0.0 


1 -t- 1 
0.1 ic 1.0 


1 7 + n Q 
1.1 m u.y 


olViiVl JloiivS^.U — Uliyoo 


Q Q 
0.0 


Q Q 
0.0 


o.D 


1 n 
i.O 


^ 1 /I Q 

< 14.8 






z. 


y. / 


^•0-1. 5 


Z.4 it 1. 1 


bMM J1311o4. 1 — 012021 


3.3 


3.2 


4.0 


1.0 


< 12.4 






2.5 


6.5 




2.6 ± 0.9 


bMM Jlcilloo. 1 — 012018 


3.3 


4.9 


4.2 


1.6 


< 17.6 






2.5 


3.8 


3.9 ± 0.5 


3.9 ±1.3 


RX JI347.5-II45 
























SMM J134728. 0-114556 


3.0 


15.5 


5.7 


3.1 


98.7 


12.8 


29.6 


2.5 


3.0 


3.0 ± 0.3 




MS1358+62 
























SMM J135957.1+623114 


3.2 


6.7 


4.4 


1.3 


< 25.5 






2.5 


1.5 


1.5 ±0.08 




A 2204 




















3 4+0.5 




SMM J163244.7+053452 


3.2 


22.2 


4.9 


5.7 


< 219.5 






2.5 


3.4 




A 221 8 
























SMM J163541. 2+661144 


2.6 


10.4 


7.5 


1.4 


53.4 


3.5 


16.0 


2.5 


1.7 


1.7 ±0.07 


9.5 ± 1.4 


SMM J163550.9+661207 


2.4 


8.7 


11.5 


1.1 


22.9 


5.9 


6.9 


2.515 


9.0 




8.4 ±0.8 


SMM J163554.2+661225 


2.3 


16.1 


21.7 


1.6 


46.4 


12.4 


13.9 


2.515 


22 




15.9 ±0.7 


bMM J 163555. 2+d6123o 


2.2 


12.8 


16.9 


1.5 


31.8 


8.3 


9.5 


2.515 


14 




12.5 ± 0.8 


SMM J163555. 2+661150 


3.3 


3.1 


3.8 


0.7 


17.1 


4.7 


5.1 


1.034 


7.6 


_ +13.0 

' -1-2.8 


2.4 ± 0.9 


SMM J163555.5+661300 


2.2 


11.3 


15.8 


1.3 


< 11.8 






2.5 


3.4 


A 9+1. u 


11.1 ±0.7 


SMM J163602.6+661255 


3.3 


2.8 


3.5 


0.6 


< 14.5 






2.5 


1.8 


1.8 ±0.08 


2.1 ±0.9 


SMM J163605.6+661259 


3.1 


5.2 


4.9 


0.9 


< 17.4 






2.5 


1.5 


1.5 ±0.03 


4.4 ± 1.1 


SMM J163606.5+661234 


3.1 


4.8 


4.6 


0.8 


< 17.4 






2.5 


1.6 


1.7 ±0.01 


4.0 ± 1.1 


A 221 9 
























SMM J164019.5+464358 


2.9 


10.0 


5.8 


2.0 


53.4 


5.0 


16.0 


2.5 


1.2 


1.2 ±0.02 




SMM J164025.5+464255* 


3.7 


5.1 


3.1 


1.7 


29.4 


2.9 


8.8 


2.5 


1.5 
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< 37.9 
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2.1 ±0.1 
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4.0 
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< 13.7 












2.0 ± 1.1 


SMM J120522.1-074431 


3.3 


3.5 


3.9 


0.9 


< 14.0 












2.7 ± 1.0 


SMM J120523.1-074516 


4.0 


3.4 


3.4 


0.9 


< 16.7 












2.3 ± 1.2 


SMM J120525.1-074512 


3.3 


4.0 


4.3 


0.9 


< 23.8 












3.3 ± 1.0 



section 13. 3|) . Taking into account that the multiply-imaged 
sources each count as one, the total number sources in the 
number counts analysis is 40. 

The first approach: For a given flux level S, only the sur- 
veyed area where 3(t < S is considered. The A^(> 5") is then 
the number of sources with > S within that area divided by 
the source plane area, Q,<s- For the cluster fields we use the 
fiuxes and sensitivity maps corrected for the gravitational 
lensing, as described in the previous section. Upper and 
lower errors are calculated using Poisson statistics weighted 
by the area, Sl<s. We use the t ables for confidence limits 
on small numbers of events from iGehrelsl l|l986D . It should 
be noted that at each S only a small number of sources is 
counted, in particular at the faint and bright end, as is re- 
flected in the error bars. The resulting number counts for 
the 850 /xm observations are plotted in Figure [S] The num- 
ber counts, A'^(> 5"), the number of sources for each data 



point and the area, n<s, are given in Table ID Due to the 
non-uniform sensitivity across the observed fields, the area 
r2<s varies with S. Consequently, A''(> S) is not uniformly 
decreasing with S, which would otherwise be expected for 
cumulative number counts determined in fields with uniform 
sensitivity. We note that even though the area of the circle 
with diameter 15" is 0.049 arcmin^, and that the counts for 
'S'850 < 0.2 are calculated from an area that appears to be 
smaller than the beam in the image plane, it should be re- 
membered that the area is in the source plane, which due to 
the large magnification would appear as a much larger area 
in the image plane corresponding to many beams. 

We use a second, alternative approach to estimate the 
number counts: 

iv(> s) = r ^ds' , (1) 
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Figure 5. The 850 /im number counts, N{> S), as determined from this work. The number counts are determined based on the data 
from 12 cluster fields and one blank field. The data from the cluster fields have first been corrected for the g ravitational lensin g. Seven 
sources have been detected with sub-mjy fluxes. This doubles the number of such sources known (compare ICowie et al.|[2002l) . Upper 
panel: The black dots show the results from the first approach for determining the number counts and the solid (magenta) line bracketed 
by the dashed show the results of the second approach (for details see text). Lower panel: The number counts from other fields have 
been plotted together with the number counts from this work. The short dashed line shows the best-fit double powerlaw, Eq. [3] and the 
long dashed line shows the best-fit Schechter function, Eq. |4l 



where n{S') is the number density of sources within the flux 
interval S' and S' + dS' . fief f is the effective area over which 
the survey is sensitive to sources with flux > 5* is given by 

^e// = fl'totat ■ C{S) . (2) 

fltotai is the total area of the survey. C{S) is the complete- 



ness function, that also takes into account the effects of the 
lensing, which is determined through simulations as follows. 
Sources are placed at random in the source plane and run 
through the lensing models (Sect.O to determine their mag- 
nification. This is done for each flux level S and compared 
with the sensitivity map to determine how large a fraction of 
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Table 2. The area in the source plane at redshift z = 2.5, an 
estimate of the source plane confusion limit (also see subsection 
14.21 1. and the area- weighted la sensitivity in the source plane. 



Cluster 


arcmin^ 


^conf 

mJy 


^850 
wghtd 

mJy 


C10016+16 


3.1 


1.7 


1.48 


A478 


3.3 


1.7 


1.65 


A496 


2.7 


1.6 


0.98 


A520 


2.2 


1.4 


0.78 


MS1054-03 


11.9 


1.8 


1.17 


A1689 


0.3 


0.4 


0.13 


RXJ1347.5- 


1145 1.2 


1.0 


1.06 


MS1358+62 


2.0 


1.3 


1.17 


A2204 


1.3 


1.1 


2.23 


A2218 


2.9 


1.2 


0.50 


A2219 


2.3 


1.4 


1.06 


A2597 


1.7 


1.3 


0.78 



^s.pl is the area of the source plane. 

Sconf is the flux confusion limit (see subsection l4.2ll . 

(T*^", , is the area-weighted 1(7 sensitivity in the source plane. 



Table 4. The 850 /xm number counts. 



5850 


N{> S) 


Nsrc ^^<3a- 


mJy 


arcmin"^ 


arcmin'^ 



0.13 


30.8+™:! 


1 


0.033 


0.20 


24.lt?5.5 


2 


0.083 


0.32 


1 fj q+16.3 


2 


0.16 


0.50 




1 


0.26 


0.80 


q 1 =-1-4.16 


2 


0.63 


1.26 


1 70-1-1.73 


3 


1.69 


2.00 




3 


4.54 


3.17 


75+0-26 
'J- '''-0.20 


14 


18.8 


5.02 


22+°-" 
"•^^-0.076 


8 


36.2 


7.96 




2 


41.2 


12.6 


024+° '"''* 


1 


42.2 



Table 5. The resulting parameters from fits to the 850 fira num- 
ber counts. 



the simulated sources would be detectable. For the sources 
from the deepest fields we use their deboosted fiuxes. Sim- 
ilarly to above, the upper and lower errors are calculated 
using Poisson statistics weighted by the area and use the 
tables fo r confi dence limits on small numbers of events from 
iGehrelsl | 1986I ). This approach allows us to take into ac- 
count the deboosted fluxes of the deepest maps. The results 
are also plotted in Figure [5] As opposed to the previous ap- 
proach, the error bars are very small at the faint fluxes due 
the larger total number of sources. We caution that these 
error bars do not contain systematic errors caused by e.g. 
uncertainties in the lensing correction. 

The number counts are probed to the faintest source in 
the survey, which has a lensing corrected flux of 0.11 mJy. 
The faint end of the number counts is dominated by the two 
cluster fields A1689 and A2218, which on the other hand do 
not contribute much at the bright end. A tentative analysis 
of the number counts show that the counts are not well- 
described by a single power-law function, but are better de- 
scribed by a double power-law 



AN{> S) 
AS 



iVo 
So 



(3) 



or another function with a turn-over, such as a Schechter 
function 



djV(> S) ^ No (S_ 
dS So [ So 



exp 



-S 

So 



(4) 



We have performed a simple x^-analysis for these two func- 
tions, as these two have been used previously for numerical 
modeling or analysis of SCUBA number counts. For this we 
have al so included the num ber counts from the SHADES 
survey l|Coppin et al.|[2"006l ) as these provide a better con- 
straint at the brighter end. We have also included the ad- 
ditional constraint that the integrated fight well below 0.1 
mJy should not be l arger the extragala ctic background fight 
jPuget et al.l Il996l : iFixsen et al] 119981 ). We note that the 
slope at the faint end is diverging and if it was to continue 



Function 


No 


So 


a f3 


Eq.|3] 
Eq.H 


658 ± 48 
1039 ± 69 


9.60tO:?° 
4.30 ± 0.08 


r,ir,-|-0.14 099 + 0.51 

^•^^-0.08 "-^^-0.34 

-2.62 ±0.10 



to much fainter fluxes it would result in an overproduction 
of the background light. The resulting parameters from this 
analysis are given in Table [5] as well as the best fit is over- 
plotted on top of the number counts shown in Figure [5] 



7.2 Fluctuation analysis 

We have performed a fiuctuation analysis, or P{D) analysis, 
on the NDF, A1689 and A2218 fields. This has previously 
been done for blank field (sub) mm data by iHughes et al.l 
(|l998h and lMalonev et all l|2005h as a statistical method to 
probe the number counts fainter than the sensitivity limit 
of the data. We measure the pixel distribution from simu- 
lated maps, which were created using an input source dis- 
tribution, convolved with the beam and added to the Monte 
Carlo maps from Section |3] As number counts for the input 
source distribution we used the Schechter function, Eq. [l] 
stepping through the three different parameters. The po- 
sitions of the simulated sources were drawn from a set of 
random positions with a normal distribution without tak- 
ing into account clustering. For the A1689 and the A2218, 
we calculated the effects of the gravitational lensing of the 
simulated sources, i.e. the magnification and position in the 
image plane. We caution that gravitational lensing will in- 
troduce similar uncertainties and effects as when calculating 
the lensing effects for the real sources and hence adding an 
extra complication for the interpretation of the results. The 
P{D) of the simulated data is compared to that of the real 
data for the three fields to determine the parameter set of 
the number counts that best fit the real data. The resulting 
parameters are listed in Table |6] along with some calculated 
counts values for comparison with the number counts deter- 
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Table 6. Results from the fluctuation analysis. 



Field 



A^o So 



N{> 2mJy) Af(> Imjy) N{> 0.5mjy) Af(> 0.2mjy) 



NDF 550 4.0 -3.75 1750 deg" 



A1689 650 4.0 -3.5 
A2218 750 6.5 -2.75 



8700 deg 



0.49arcmin ^ 2.4arcmin 



24600 dcg-2 114600 deg-2 

6.8arcmin~^ 31.8 arcmin"^ 

30700 deg-2 72400 deg-2 

8.5arcmin~^ 20.1 arcmin"^ 



mined in the previous subsection. As seen for both the NDF 
and A1689, a < —3.5, gives a very steep function especially 
at the faint end which would strongly overproduce the EBL. 
The resulting parameters from A2218 are reasonably close 
to those deduced the Schechter function fit to the number 
counts. We note that although all fields are roughly equally 
deep in the image plane, the faint counts are best probed 
by A2218, since NDF is not gravitationally lensed so that 
that field does not probe faint fluxes very well, and A1689 
covers only a very small area in the source plane. As P{D) 
analysis is a statisti cal tool it is be s t app lied on large fields 
as was done by e.g. iMalonev et al.l (|2005l ). 



7.3 Comparison with other surveys 

Here we will compare the derived number counts with those 
determined through other surveys, both lensed surveys and 
blank field surveys. The number counts from other studies 
have been plotted together with the number counts from 
this work in Figure [5] 

Lensing surveys: Three other studies of SCUBA ob- 
servations of cluster fields have been published: For the 
UK-SCUBA Lens Survey seven cluster fields were targeted 
and number counts were determined to Ss50 = 0.5 mj y 
(e.g..lSmail et al.lll997l : iBlain et al.ll 19991 : ISmail et al.ll2002h . 
ICowie et al.l ( 20021 ) obtained deeper SCUBA observations 



of three of t hese fields. A shallowe r cluster survey was per- 
formed by (|Chapman et al.ll2002l ). in which eight clusters 
were observed with SCUBA, however, no sub-mjy sources 
were dete cted. For > 1 mJy, our number counts agree with 
those of (jChapman et al.l |2002| ). Here we will focus the 
comparison on th e surveys from ICowie et al.l (|2002| ) and 
ISmail et all l|2002l ). As both those surveys are relatively 
small in area, a comparison is only interesting where such 
surveys hav e thei r strength, namely at the faint fiuxes. 
ICowie et al.l l|2002l ) detect five sub-mJy sources. We detect 
seven sub-mJy sources, and thereby double the number of 
known sub-mJy sources. Our faint number counts are in 
good agreement with those of Cowie et al. and Small et al. 

Blank field surveys: Blank field surveys are surveys with 
no strongly gravitationally lensing clusters present in the 
surveyed area. Such surveys typically cover much large ar- 
eas than the lensed surveys, and are limited in depth by the 
blank field confusion limit (~2.0mJy). Hence the strength 
of those surveys lies at brighter fluxes. Several such surveys 
have been carr ied out: The Canada-UK Deep SCU BA Sur- 
vey (CUDSS) (|Eales et al.ll2000l : IWebb et al.ll2003l ). covered 
75arcmin^ to the blank field confusion limit. The 8m Jy- 
survey (Scott et al.) covered an area of 260 arcmin^ to a 



flux limit >5mJy/beam. The Hubble Deep Field North 
(HDF-N) has been surveyed extensively, which has been 
brought to gether in the so-c alled "HDF-N SCUBA su- 
permap bv lBorvs et ai] (|2003l ) , which covers 165 arcmin'^ to 
depths between 0.4 and 6 mjy/beam. iBarger et al.l (|l999l ) 
have surveyed the Hawaii Survey Fields covering an area 
of 104arcmin^ to a flux limit of 8mJy with a small area 
of 7.7 arcmin^ almost to the confusion limit. ICoppin et al.l 
(|2005l ) surveyed 70 arcmin of the Groth Strip to a depth of 
Icr rms. « 3.5 mJy. A re -analysis of several blank field sur- 
veys was carried out bv lScott et al.l ([2006). Fi nally the re- 
sults of the SHADES survey were published bv lCoppin et al.l 
(,2006). where 720 arcmin^ were covered to a noise level of 
~ 2 mJy uncovering more than a 100 submill i metre galax- 
ies. Through a detailed analysis ICoppin et al.l (|2006l ) deter- 
mine the first differential submm number counts and fit their 
results with a double power law. With minor deviations, 
there is an overall good agreement between the bright end 
the number counts of the work presented here and previous 
work. Though, we do find that the slope of the power-law 
at the bright end is a bit steeper than previous work (a ~ 
1.9-2.2). 



7.4 Resolving the extragalactic submillimetre 
background light 

Using the differential number counts from equation O we 
calculate the integrated background light. At S ~ 0.10 mJy, 
the integrated background produced by our sources is 
comparable to the backgrou nd light detected with COBE 
(|Puget et al.lll996l : [Fixsen et al.lll998t ). Given that sources 
with fluxes below 0.1 mJy also contribute to the integrated 
background, there is a possibility that our counts overpredict 
the integrated background somewhat. The overproduction of 
the background light, which is caused by the shape of the 
number counts at the faint end, is possibly due to the low 
number statistics. Given the differential counts from equa- 
tion [3] the dominant contribution to the integrated back- 
ground light comes from the sources with fluxes Sgso be- 
tween 0.4 mJy and 2.5 mJy with 50 per cent of the back- 
ground resolved at 1 mJy. The lat t er is in agreement with 
the results frornlSmail et all l|2002l ). ICowie et all (|2002l ') and 
IChapman et al.1 (|2002D ~ Sources with Sgso > 2.5 mJy con- 
tribute ~ 25 per cent to the integrated background, of which 
sources with ^gso > So contribute only ~ 10 per cent. This 
means that the bulk of the submm energy output from the 
submm galaxy population arises from sources just fainter 
than the blank field confusion limit. 
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8 CONCLUSIONS 

We have conducted a deep submm survey using SCUBA. We 
have observed twelve clusters of galaxies and the NTT Deep 
Field. The total area surveyed is 71.5 arcmin'^ in the image 
plane. For the cluster fields the total area in the source plane 
is 35 arcmin^. This is the largest deep submm lens survey of 
its type to date. The gravitational lensing reduces the 
confusion limit allowing for observations of sources with 
Ssso < 2 mjy. The data have been analysed using Monte 
Carlo simulations to quantify the noise properties, Mexican 
Hat Wavelets have been used for source extraction and simu- 
lations were performed to quantify the error of the analysis. 

• In total 59 sources have been detected, of which 10 have 
flux densities below the blank field confusion limit. Four are 
associated to cluster cD galaxies. Three sources in the field of 
A2218 are multiple images of the same galaxy, and in A1689 
five are associated with two multiply-imaged galaxies. 

• The number of sub-mjy sources is seven, which doubles 
the number of such sources. 

• The integrated number counts are probed over two 
decades in flux down to 0.10 mJy. The number counts cannot 
be described by a single power law, but have to be described 
by a double power law or another function with a turn-over. 
Describing the differential counts by a double power law 
function we find that the turn-over is ~ 6 mJy. At 1 mJy 
the number counts are ~ 10* deg~^ and at 0.5 mJy they are 
^2x10* deg~^, based on derived differential counts. 

• Another key result is that essentially all of the inte- 
grated submm background is resolved. At 1 mJy 50 per cent 
of the background is resolved, and at 0.4 mJy 75 per cent 
is resolved. The dominant contribution to the background 
comes from sources with fiuxes Ssbo between 0.4 mJy and 
2.5 mJy, while the bright sources with fluxes Sgso > 6mJy 
contribute only 10 per cent. This means that the bulk of the 
energy comes from submm galaxies with fluxes just below 
the blank fleld confusion limit. 

The submm number count distribution is an observ- 
able for the submm galaxy population as a whole and pro- 
vides strong constraints on the models describing the submm 
galaxy population and their evolution. While the submm 
number counts are well studied at the bright end (> 2mJy), 
the faint end (< 2mJy) and the extremely bright end 
(> 20mJy) remain difficult to probe. The extremely bright 
end is challenged by the steep counts. The faint end is chal- 
lenged by the blank fleld confusion limit. The present survey 
has made a substantial contribution to the faint end, how- 
ever, it essential to follow this through with future instru- 
mentation such as more sensitive instruments like SCUBA-2 
and LABOCA, with which an even larger number number of 
strongly lensing clusters can be surveyed, larger telescopes 
such as the LMT and CCAT for which the blank field con- 
fusion limit will be lower, and of course ALMA which will 
be able to study the faint sources in great detail. 
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APPENDIX A: DESCRIPTION OF THE 
INDIVIDUAL FIELDS 

C10016+16 In this field, four point sources were detected 
with 3 < S/N < 4. The magnification factors for the 
four sources are between 1.2 and 2.2. The corrected fiuxes 
of the sources are between 2.4 and 5.2 mJy. We note 
that three sources have more than 5 per cent of their 
posterior probability distribution below mJy. None of the 
sources have significant 450 /im fiux detecti ons. A shallower 
SCUB A map of this field is presented in IChapman et al.l 
(|2002| ). where the depths is about a factor two shallower 
than our map. Within the positional uncertainties the 
two sources from Chapman et al. are coincident with 
SMMJ001834.2+162517 and SMM J001835. 1+162559, 
though we find that the observed fluxes are fainter than in 
Chapman et al. The mass model is based on the results from 
Natarajan (private communication) with a substantially 
more detailed description compared to that of Chapman et 
al. 

A^^S T his cluster i s well-known in cooling-flow stud- 
ies (e.g., IWhite et af] [1991). Some of the SCUBA data 
presented here were obtained by others to study the cooling- 
flow, however, the cooling-flow has not been detected in the 
data. Four point sources are detected. With a detected flux 
of Ssso = 25±3mJy the source SMM J041327.2-H102743 is 
the brightest source in the survey. This source has been 
studied in detail and is identifled with a type-1 quasar 
at redshift z — 2.837 (|Knudsen et al. I l2003l , there denoted 
SMMJ04135-F-10277). The three fainter sources have S/N 
between 3.3 and 4.2. All four sources have magnification 
factors of 1.2-1.3. The fluxes of the three other sources are 
5.6-7.3 mJy. We note that two sources have more than 5 
per cent of their posterior probability distribution below 
mJy. The quasar is the only source in this fleld with 
450 i-im flux detection. Close to the SE-edge a fifth bright 
source is detected, however, as it is less than 1.5 beam 
from the edge it is not included in the catalogue. The 
mass model is a simple model, which we constructed based 
on the published velocity dispersions: the model includes 
the cluster potential and the potentia l of the cD galaxy 
(|Zabludoff et al.lll990l : lAllen et al.|[l993D . 

A496 This is the lowest redshift cluster in the sur- 
vey. Nine point sources have been detected. Even though 
this cluster has not been observed to the blank field con- 
fusion limit, the large number of sources might introduce 
extra uncertainties on the derived parameters. Three 
sources towards the centre of the fleld are just 14" (just 



14 K.K. Knudsen, P.P. van der Werf, J. -P. Kneib 




R.A, (3000.0) 

Figure Al. The signal-to-noise (S/N) 850 /^m SCUBA maps of the clusters C10016+16, A478, A496 and A520. The overlayed contours 
represent S/N = 3,4,5,6,7, but for A478 the contours represent S/N = 3,4,5,6,8,10,12,14. The stars indicate the position of the detected 
sources. 



smaller than a beam) from one another. One of the central 
sources, SMM J043337.8-131541, is coincident with the 
cD. The two other central sources, SMM J043337.4- 131558 
and SMM J043336. 5— 131547, are so close (just less than a 
beam) to the centre of the cD galaxy, that they are likely 
associated with cD galaxy. The latter of those two sources 
has a probable detection of 450 ^m emission. All three 
central sources will be excluded from the further analysis in 
this paper. In Figure IXT] source SMM J043338.9- 131444 
has no S/N = 3 contour as it is located in a depression in 
the background. The six sources not associated with the 
central cD galaxy are magnified by factors 1.3-1.4 and have 
unlensed fluxes of 3-6.3 mjy. We note that two sources have 
more than 5 per cent of their posterior probability distribu- 
tion below mJy. This large number of sources in a low-z 
cluster field is surprising. Follow-up observations indicate 
that they are not cluster members, which might otherwise 
be expected because of the low redshift of the cluster. Like 
A478, the lens mass model is a simple model including the 
cluster potential and the cD galaxy (|Peletier et al.l Il990l : 
[Zabludoffct al. IQqJ). 

A 5 20 The optical centre and the X-ray centre of A520 
is not coincident, and the cluster seems to be undergoing 
strong dynamical evolution, as the cD galaxy is not l ocate d 
at the centre of the X-ray emission (|Proust et al]|2000l ). 
Our SCUBA map is about 1' E of the X-ray centre 
ijCovoni et al.ll200ll l. Four point sources have been detected 



with lensing corrected fluxes between 0.7 and 3.4 mJy. 
We note that two sources have more than 5 per cent of 
their posterior probability distribution below mJy. The 
brightest source in the field, SMM J045409. 7-1-025510, has 
a possible detection of 450 fim flux. T he mass model is 
based on the genera l cluster potential (jWhite et al.l Il997l : 
ICarlberg et alii 19961 ). 

MS1054—03 For this cluster the deepest multi- wavelength 
data set exists, ranging from radio to X-ray. It is a part 
of the FIRES-pr oj ect (Faint IR Extragalactic Survey; 
iFranx et al] I2OO0I : iForster Schreiber et all 120061 ). which 
includes the deepest near-IR imaging of a cluster taken 
with ISAAC at VLT. The area covered with ISAAC and 
other instruments is about 5'x5', so we decided to obtain 
three pointings to cover a larger area of the field and take 
advantage of the excellent data available for follow-up 
studies. The three pointings, which are denoted by 5", 
A'^ and NW , according to the relative position cover a 
signification part of the FIRES field. The S pointing is 
centered at the cluster centre. We detect nine sources with 
fluxes between 3.5 and 5.0 mJy. We note that one source 
has more than 5 per cent of their posterior probability 
distribution below mJy, however, this source has already 
been identi fied with a Dist ant Red Galaxy at redshift 
z = 2.423 (|van Dokkum et a l. 2004; Knudsen et"all |2005| ) 
and thus we do not consider this a spurious detection. 
The largest fraction of the sources are located in the N 
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Figure A2. The signal-to-noise 850 /im SCUBA map of the cluster MS1054-03. The overlayed contours represent S/N = 3,4,5,6,7,8. 
The black boxes indicate the positions of the detected sources. To cover a large area of the existing multi- wavelength data (see text for 
details), we have obtained three pointings for this field. 



pointing, while the NW pointing is a lot sparser. This 
suggests a level of clustering of SCUBA sources, though 
the field is too small for a reliable clustering analysis. Our 
map is about three times deeper than the shallo w map of 
the a rea of the ^-pointing published by Chapman et al.l 
||2002D . They find one source, which is off-set ~ 25" north 
of SMMJ105703. 7-033730 and not at all detected in 
our much deeper map. The gravitational lensing is not 
particularly strong for this cluster, which is partly related 
to the relative high redshift. T he mass model i s based 
on the overall cluster potential (|Tran et al.l 1 19991 : P. van 
Dokkum, private communication). 

A 1689 With ~ 34 hours raw integration time for a 
single pointing and a la r.m.s. ~0.7mJy, this is one of the 
deepest maps of the survey. A1689 is a cluste r known to 
have an exceptionally large Einstein radius fe.g.. lKing et al.l 
|2002| ). Because the gravitational lensing is so strong, the 
source plane area at redshift z = 2.5 is only 0.3arcmin'^, 
i.e. 20 times smaller than the image plane area or the 
field of view of SCUBA. We detect nine SCUBA sources, 
and note that this field might be suffering from confusion 
due to the large number of sources. The three central 
sources are approximately one beam from one another, and 
the same is the case for two eastern sources. The sources 
have observed fluxes between 2.6 and 5.4 mJy. When 
correcting for the gravitational lensing magniflcation the 
fluxes are between 0.11 and 1.3 mJy. The central source 
SMMJ131129.1-012049 hasa probable detection of 450 
flux. 



Two multiply-imaged galaxies have been identified among 
the submm galaxies in this field. SMM J131129. 1-012049 
and SMM J131134.1-012021 have been identified 
with the triple -imaged system 5 (as numbered in 
iBroadhurst et al.l (j2005)), while SMM J131132. 0-011955, 
SMMJ131129.8-012037 and a small contribution to 
SMMJ131134.1-012021 arise from either system 24 or 
29, both quintuple-imaged galaxies. Recent optical spec- 
troscopy has shown that all these galaxies have redshifts 
~ 2.5 (Knudsen et al., in prep., and Richard et al., in prep.). 
A detailed analysis of the identification will be discussed 
in a future paper along with additional multi-wavelength 
follow-up. The mass model will be presented in detail in 
Richard et al. [in prep.) and Limousin et al. {submitted). 

RX ,11347.5— 1145 This is the most X-ray-luminous 

cluster known ( jAllen et al.l I2OO2I ) . In this field we detect 
one source, SMM J134728.0-114556, which has Ssao = 
16.2±3.1mJy and S'450 = 98.7±29.6 mJy. The source is 
strongly lensed and has an unlensed flux of 4.5 mJy. 
Furthermore, a large, extended source near the cluster 
centre is present. This is the Sunyaev-Zel'dovich effect 
re ported in Komat su et al.l l| 19991 ). The mass model is based 
on lCohen fc Kneitr(|20oi v 

MS1358+ 62 This map is relatively empty with only 
one detection SMM J135957.H-623114, with S'sso = 
6.8±1.3mJy; 4.4 mJy after correcting for the lensing. 
These data were first obtained to study the str ongly lensed, 
redshift z = 4.92 galaxy MS1358-I-62-G1 (|Franx et all 
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Figure A3. The signal-to-noise (S/N) 850 /^m SCUBA maps of the clusters A1689, MS1358-I-62, RX J1347.5-1145 and A2204. The 
overlayed contours represent S/N = 3,4,5,6,7,8. The stars * indicate the positions of the detected sources. 



1997 ) , which however was not detected (|van der Werf et al.l 
200 ll ). We here find a 3a upper limit for Gl of Sgso < 



66°14 



4.8 mjy. A detailed mass model describes the po tential for 
this cluster (|Franx et al.lll997l : [Saiitos et al.ll2004l '). 



the shallowest SCUBA ob- 
survey. One point source. 



A 2204 This field has 
servations of the whole 
SMMJ163244.7+053452, has been detected at a S/N = 4.9 
with an observed fiux of Ssso = 22.2±5.7mJy making it 
the second brightest source in the catalogue. This source is 
lensed by more than a factor 3, resulting in a corrected fiux 
of ~ 7 m Jy. 

A 22 18 Together with A1689 and the NTT Deep Field, this 
field is the deepest data taken for the survey. The data for 
this field cover an area corresponding to more than two 
pointings. The data was used as a case study for the source 
extraction method, the Mexican Hat Wave lets algorithm, 
applied for this survey (|Knudsen et al.|[2006l ). 

In this field nine sources were detected. The three 
source, SMM J163550.9+661207, SMM J163554.2-)-661225 
and SMM J163555. 2-1-661238, have been identified as 
the same, multip ly-image source at redshift z — 2.516 
l|Kneib et al.ll2 Q04h. and is referred to as SMM J16359-h6612. 
The source SMM J163555. 2-1-661150, which is detected both 
at 850 jxm and 450 /xm, is coincid ent with a known galaxy, 
#289, with redshift z = 1.034 (|Pello et al l 1 19921) and is 
also detected at 15 /xm with ISOCAM (e.g., Metcalfe et al.l 
l2003h . The relatively bright source SMM J163541. 2-^661144 
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Figure A4. The signal-to-noise 850 SCUBA map of 
the cluster A2218. The overlayed contours represent S/N = 
3,4,5,6,8,12,18,24. The black boxes indicate the positions of the 
detected sources. 



is detected at both 850 /im and 450 fim. The sources in 
A2218 have observed fiuxes between 2.8 and 16.1 mJy. 
The lensing correc ted fluxes a r e 0.4 - 6.1 mJy. The mas s 
model is ba s ed on iKneib et al.1 (|l996l l: lEUis et all (|200lf ): 
iKneib et all |2004h . 
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Figure A5. The signal-to-noise (S/N) 850 /xm SCUBA maps of the clusters A2219 and A2597 and of the blank field NTT Deep Field. 
The overlayed contours represent S/N = 3,4,5,6. The stars * indicate the positions of the detected sources. 



A 22 19 In this field we detect two point sources that 
both have possible detec tions of 450 /xm flux. T his field 
was also a part of the IChapman et al.1 (j2002h sample, 
though their observations are shallower. The source 
SMM 164019.5-1- 4 64358 agrees well with their finding. In 
IChapman eraP (|2002[ l the source C SMM J16404-f 4643 
has an upper limit at 850 fim, while being detected at 
450 /xm. We do not get a significant detection of the source, 
but MHW does suggest a S/N ~ 1.9 detection, which 
corresponds to a fiux of ~ 2.3 mjy; the 3cr upper limit is 
< 3.6 mJy. We note that one source has more than 5 per 
cent of its posterior probability distribution below mJy. 
For the source B SMM J16403+46437 MHW suggests a 
S/N ~ 2.3 detection. The source D SMM J16404-f 4644 is 
within the edge region which have trimmed from the map. 
In our map a positive fluctuation is present, though it 
does not have the characteristics of a significant 850 ^m 
detection. Furthermore, MHW also finds a S/N ~ 2.9 
source at a, 5^JJ^W^2^+46°42'25". The mass model is 
described in lSmith et"aLl l|2005h . 

A2597 In this field, two point sources were detected with 
S/N > 3. The brightest source, SMM J232519.8-120727, 
is a 12 mJy source located in the centre of the map and 
is coincident with the c D galaxy of the c l uster , which is 
a well-known AGN (e.g., iMcNamara et al.lll999l l. The cD 
galaxy is excluded from the rest of the analysis in this 
paper. The other source, SMM 232523.4-120745, is a 5 
mJy source, which also has detected 450 fj,m flux. We note 



that this source has more than 5 per cent of their posterior 
probability distribution below mJy. Hence, if this is indeed 
a spurious source, then no high-z background sources were 
detected in this field. The mass model includes both the 
overall potential of the cluster and that of the cD galaxy 
(jSmith et al.lll990l : IWu fc Xuell2000l ). 

A^rr Deep Field Th is field, the blank field of the 
survey (Ar nouts et al.1 [l999) . is one of the deepest fields 
of the survey. Five sources have been detected with fiuxes 
between 3 and 4 mJy. None of them have detected 450 /xm 
fiux. In Figure HI source SMM J120520.6-074448 has no 
S/N = 3 contour as it is located in a depression in the back- 
ground. A 1.2 millimeter map of the NTT Deep Field has 
been obtained with the Max-Planck Millimetre Bolometer 
(MAMBO) covering a larger area of t he NTT Deep Field 
than the S CUBA map presented here jPannerbauer et al.l 
|2002| . [ 20041 ). Considering that 1.2 mm probes a part of the 
modified blackbody where the flux is fainter compared to 
the 850 /xm, the MAMBO map is a bit shallower than the 
deep SCUBA map. Two MAMBO sources are covered by 
the SCUBA map. The source MM J120517-0743.1 is very 
close to the edge of the SCUBA map, where there are no 
indications of a source. The source MM J120522-0745.1, 
which has a radio detection, is only 6" from the submm 
source SMM J120523.1-074516. The radio detection is 
coincident with the submm source. 
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